YITP-SB-06-52 

UB-ECM-PF-06-38 

MADPH-06-1474 



Gamma Ray Burst Neutrinos Probing Quantum Gravity 

M.C. Gonzalez-GarciaQ 
Institucio Catalana de Recerca i Estudis Avangats (ICREA) 
Departament d'Estructura i Constituents de la Materia 
Universitat de Barcelona, 
647 Diagonal, 08028 Barcelona, Spain 
and 

C.N. Yang Institute for Theoretical Physics, 
SUNY at Stony Brook, 
Stony Brook, NY 11794-3840, USA 

Francis HalzerQ 

Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

Abstract 

Very high energy, short wavelength, neutrinos may interact with the space-time foam predicted 
by theories of quantum gravity. They would propagate like light through a crystal lattice and be 
delayed, with the delay depending on the energy. This will appear to the observer as a violation of 
Lorenz invariance. Back of the envelope calculations imply that observations of neutrinos produced 
by gamma ray bursts may reach Planck-scale sensitivity. We revisit the problem considering two 
essential complications: the imprecise timing of the neutrinos associated with their poorly under- 
stood production mechanism in the source and the indirect nature of their energy measurement 
made by high energy neutrino telescopes. 



* Electronic address: lcoiicha@insti. physi cs. sunysb.edu| 
Electronic address: .flhalzen@facstaff.wisc.edul 



1 



I. INTRODUCTION 



It has been realized for some time that particle physics experiments using cosmic beams 
may be sensitive to Planck scale physics[l|]. Neutrino telescopes in particular have un- 
matched sensitivity to violations of Lorenz invariance (LIVs) and the equivalence principle 
by exploiting the high statistics observations of the guaranteed high energy atmospheric 
neutrino beam in conjunctions with methods demonstrated by the Superkamiokande and 
Macro collaborations [2|]. LIVs are the sources of observable flavor oscillations caused by the 
different speed of light of the three neutrino flavors. While the atmospheric neutrino data 
will reach PeV energy, the observation of cosmic beams may extend the neutrino sample to 
EeV energy. Of particular interest in this context is the observation of high energy neutrinos 
from gamma ray bursts (GRBs) whose short-time burst nature make more direct tests of 
Lorenz invariance possible. Their energy is such that the measurements may be sensitive to 
modifications of Lorenz symmetry associated with Planck scale physics jsl, Q, 

Lorenz invariance can be tested by measuring the travel times of particles over cosmic 
distances. The method is familiar from pulsar observation [gI. The pulsar provides a clock 
at the time of emission. It has recently been suggested that more sensitive observations can 
be made by measuring the relative time of photon and very high energy neutrino emission 
by GRBs The photon emission in a burst lasting seconds starts the clock for a neutrino 
whose early or delayed arrival time At will be a signature of LIV. Note that at the energies 
relevant to this problem neutrinos are, like photons, massless particles. 

The measurement conceptually tests, in an essentially model-independent way, LIV as- 
sociated with the possible breakdown of Lorenz symmetry at the Planck scale. In some 
theories unifying quantum mechanics and gravity the continuous spacetime of general rel- 
ativity becomes quantized in volumes that cannot be further subdivided. Particles with 
very high energy and small wavelength may sense this spacetime foam and be delayed, just 
like photons traveling through a crystal. The phenomenology can be discussed as a generic 
modification of the dispersion relation 

E'=p' + m'±E\—^r±... (1) 

Planck 

for particles with energy E « ^Epianck- The time delay of a particle with energy E emitted 
at redshift z (distance d) is 



^^E Planck Hq Jq A/fim(l + z'Y + ^ C ^Epianck 

(2) 

A sensitive test requires good timing At, high energy E and large distances d. Back-of-the- 
envelope is sufficient to conclude that GRBs provide the opportunity to reach the Planck 
scale, i.e ^ ~ 1, at least for n=l, although there is little sensitivity for higher values of n. It 
is also misleading. 
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As will become clear as the arguments in this paper develop, it will be challenging to ever 
argue that the delay of a high energy neutrino relative to the light of a GRB requires LIV 
rather than a straightforward astrophysics explanation associated with the poorly under- 
stood nature of the source. We will argue indeed that the fireball phenomenology of GRBs 
can accommodate the emission of high energy neutrinos over timescales from seconds to days 
and beyond. Emission over such timescales is actually expected, as we will argue further on. 
In particular, we will give an example in the next section how delays of very high energy 
neutrinos by more than 1,000 seconds from the start of the photon display, can be accom- 
modated by fireball phenomenology rather than by LIV. This renders the implementation 
of the proposal in reference 0] more challenging. 

We here also revisit an alternative proposal to exploit the emission of two neutrinos from 
a smg le GRB 0, Even though the neutrino rate from GRBs is expected to be order 10 
per kilometer squared per yearj^], events with two neutrinos are not rarej^. This is because 
neutrinos are mostly produced by relatively nearby and relatively energetic bursts. As an 
example we refer to the calculation of reference [9|]. The large burst-to-burst fluctuations can 
be seen from their figure 2. Although their conservative estimates predict less than 10 events 
per year, applying Poisson statistics we anticipate that bursts with 2 neutrinos should be 
observed within 2 year intervals. 

We will emphasize how GRB fireball phenomenology anticipates a very attractive timing 
of the neutrino emission, with TeV neutrinos produced before, PeV neutrinos in coincidence 
with, and EeV neutrinos after the photon display. I.e. higher energy particles come later 
and a reversal of this order can provide us with evidence for LIV. 

We finally confront the additional problem that kilometer-scale neutrino telescopes do not 
measure neutrino energy directly. For the dominant z/^ signal the detector only measures the 
energy-loss of the secondary muon when it passes through the detector [l^]. We will derive 
the probability distribution of neutrino energies that accommodates the observed energy in 
the detector. Our results on energy measurement presented should be of interest beyond 
their application to the problem considered here. The problem is less severe for electron 
and tau neutrinos where the detectors are total absorption calorimeters for the secondary 
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showers produced by the neutrino interactions 

In the end we conclude that revealing Planck-scale physics with IceCube may be possible, 
but will be challenging. Reading the time of the astrophysical GRB clock is not straightfor- 
ward and the energy of the neutrinos is not directly measured. We point out scenarios where 
the first problem can be overcome, in the meantime we will just have to wait for observations 
that will clearly define the dilemma, if any. The second problem is interesting and impacts 
all observations with high energy neutrino telescopes. A straightforward solution is to do the 
science with electron and tau neutrinos only. The energy of their secondary showers can be 
measured with a linear energy resolution of 20%. Though GRBs represent a background- free 
signal in neutrino telescopes because of the clustering in time and in direction, the event 
rates are still uncertain. They will exceed 10 events per kilometer square year, in the case 
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where the baryon loading in GRB fireballs accommodates the observed cosmic rays. With 
the low event rates anticipated doing the Planck science with showers may not be a practical 
solution and we may be forced to explore the long range of the secondary muons in order 
to collect the enhanced rates expected for muon neutrinos. Their better pointing will also 
facilitate background rejection. How well one reconstructs muon energy of individual events 
has to the best of our knowledge not been systematically analyzed. Our results should be 
of interest beyond their application to the problem considered here. 



II. NEUTRINO TIMING BY GAMMA RAY BURSTS 



We here review how GRBs produce a sequence of TeV, PeV and EeV neutrinos separated 
in time. They provide a clock with later time implying higher energy. We propose to exploit 
this feature for measuring particle delays signaling LIV. 

Gamma ray bursts are perhaps the best motivated sources of high-energy neutrinos 

m 
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131 ]. The collapse of massive stars to a black hole has emerged as the likely origin of the 
"long" GRBs with durations of 10 seconds on average. The GRB rate is consistent with the 
rate of supernovae with progenitor masses exceeding several solar masses. In the collapse 
of the massive star a fireball is produced which expands with a highly relativistic velocity 
powered by radiation pressure. The fireball eventually runs into the stellar material that is 
still accreting onto the black hole. If it successfully punctures through this stellar envelope 
the fireball emerges to produce the GRB display. While the energy transferred to highly 
relativistic electrons is thus observed in the form of radiation, it is a matter of speculation 
how much energy is transferred to protons. 

The phenomenology that successfully describes GRB observations is that of a fireball 
expanding with highly relativistic velocity, powered by radiation pressure. The observer 
detects boosted energies emitted over contracted times; without this a description of the 
extreme observations is impossible. The millisecond variations of the GRB flux require an 
original event where a large amount of energy is released in a very compact volume of order 
100 km. The dynamics of the fireball is actually reminiscent of the physics of the early 
expanding universe. Initially, there is a radiation dominated soup of leptons and photons 
and few baryons. It is hot enough to freely produce electron-positron pairs. With an optical 
depth of order ~ 10^^, photons are trapped in the fireball. It cannot radiate. This causes the 
highly relativistic expansion of the fireball powered by radiation pressure. The fireball will 
expand with increasing velocity until it becomes transparent and the radiation is released in 
the display of the GRB. By this time, the expansion velocity has reached highly relativistic 
values of order 73OO. 

With a relativistic expansion of 73OO it may be possible that protons are accelerated to 
energies above 10^° eV in the fireball that accommodates the GRB observations Esl. 16|. 
GRBs within a radius of 50-100 Mpc over which protons can propagate in the microwave 
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background, may therefore be the sources of the ultra high-energy cosmic rays. The as- 
sumption that GRBs are the sources of the highest energy cosmic rays does determine the 
energy of the fireball baryons. Accommodating the observed cosmic ray spectrum of extra- 
galactic cosmic rays requires roughly equal efficiency for conversion of fireball energy into 
the kinetic energy of protons and electrons. In this scenario the production of neutrinos of 
100-1000 TeV energy in the GRB fireball is a robust prediction because pions, and there- 
fore neutrinos, are inevitably produced in interactions of accelerated protons with fireball 
photons 0, HHQ. 



There are several other opportunities for neutrino production following the chain of events 
in the production of a GRB. 

• TeV neutrino production in the stellar envelope ~ 10 seconds before the GRB display. 

The core collapse of massive stars has emerged as the likely origin of the "long" GRB 
with durations of tens of seconds. The fireball produced is likely to be beamed in 
jets along the rotation axis of the collapsed object. The mechanism is familiar from 
observations of jets associated with the central black hole in active galaxies. The jets 
eventually run into the stellar material that is still accreting onto the black hole. If the 
jets successfully puncture through this stellar envelope they will emerge to produce a 
GRB. While the fireball penetrates the remnant of the star, the fast particles in the tail 
will catch up with the slow particles in the leading edge and collide providing another 
opportunity for pion production yielding neutrinos of tens of TeV energy. Bursts within 
a few hundred megaparsecs (~ 10 bursts per year as well as an additional unknown 
number of "invisible" bursts from failed GRBs) may actually produce large rates of 
TeV neutrino events in a kilometer-scale detector and, possibly, observable rates in a 
first-generation detector such as AMANDA[19!]. Interestingly, failed jets that do not 
emerge will not produce a visible GRB but do produce observable neutrinos [itI. fisl . 

• PeV neutrino production inside the expanding fireball for a duration of 10 seconds 
coincident with the GRB. 

The production of PeV neutrinos in cosmic ray producing GRB fireballs is a robust 



predict ion [111]. Neutrinos are inevitably produced in interactions of accelerated pro- 
tons with fireball photons. As the kinetic energy in fireball protons increases with 
expansion, a fraction of their energy is converted into pions once the protons are 
accelerated above threshold for pion production. 

EeV production in the interstellar medium over minutes days after the GRB display. 

Afterglow observations show that external shocks are produced when the GRB runs 
into the interstellar medium. The hadronic component in the fireball can be rapidly 
depleted by means of photopion process on time scales of 10^ ~ 10^ seconds. It has 
been argued that this "aftershock" explains the accumulating evidence from recent 
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SWIFT observations for a second high energy component in the GRB as well as for 



the sharp decline of the X-ray light curves that take place over similar timescales 2G ] . 
This is the phase in which the highest energies are reached, 100 EeV for cosmic rays 
and 1 EeV for the secondary neutrinos. The neutrino rates are predicted to be low in 
IceCube but represent an opportunity for radio or acoustic detectors such as RICE 
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with high threshold, but also larger effective area 

It is this last phase of neutrino production that will interfere with the proposal of 
Jacob and Piranj^. Optimistically assuming that one will make multiple observations 
and that, contrary to evidence, GRBs are somehow standard candles, one could still 
try to distinguish LIV from astrophysics induced delays by establishing the linear 
dependence of the time delay on the distance; see Eq. ([2]). 

While the arrival time of neutrinos is measured with microsecond accuracy, determination 
of the energy is a considerable challenge. Both are required to test LIV. 



III. HIGH ENERGY NEUTRINO TELESCOPES AND ly^ ENERGY MEASURE- 
MENTS 

In neutrino telescopes, neutrinos are detected through the observation of Cerenkov light 
emitted by charged particles produced in neutrino interactions. The neutrino induced events 
can be categorized as either muon tracks or showers. Cosmic ray muons and muons from 
charged current (CC) t'^ interactions are the origin of tracks. Showers results from neutrino 
interactions — Ue or ly^ CC interactions, and neutral current (NC) interactions initiated by all 
three flavors — inside or near the detector. Because of the large range of muons, kilometers 
to tens of kilometers for the energies considered here, the effective volume of the detector 
for muon neutrinos is significantly larger than the instrumented volume. Furthermore, the 
angular resolution for muon tracks is superior to that for showers. For example, for Icecube 



the angular resolution for muon tracks ~ 0.7° 22| which allows a search window of solid 
angle A^ioxio ~ 3 x 10^"^ SD, while for showers the angular resolution is significantly worse. 

While the large effective volume and the better background rejection of atmospheric back- 
grounds because of the superior angular resolution favors the observation of muon neutrinos, 
measurement of their energy is indirect. Neutrino telescopes only measure the energy E"^*^^ 
deposited by the muon when traversing the detector which is only indirectly related to the 
original neutrino energy. Particle physics determines the probability distribu- 
tion for the parent neutrino energy of a muon track that deposits an energy E'^'^^ inside the 
detector: 

^^^^1^?^) -^f'^ C '^'^ C '-k^'^^ ^» ^^^2' ' 

(3) 
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is the incident neutrino spectrum, ^^^{E,y, E^) is the differential neutrino CC interac- 
tion cross section and F{E'^, E^'^, I) represents the probabihty that a muon produced with 
energy E^ arrives at the detector with energy E^^ after travehng a distance /. / is obtained 
by propagating the muons through the ice taking into account energy losses from ioniza- 



tion, bremsstrahlung, e^e pair production and nuclear interactions 23|]. The details of 



the detector are encoded in the effective area A^jf and /^jn is the minimum muon track 
length required for the event to be detected. Norm is the normalization factor obtained 
from J^a.,f{E,\Ef^^P)dE, = 1. 

We show in FigiHthe probability distribution function f{Eiy\E^'^''^) ( for convenience we 
plot it as f (\ogiQ{E^)\E'^'^^) = \n{10) E^f{E^\E'^'^^)) as a function of the neutrino energy for 
different values of E'^'^^ and for two different neutrino spectra, a broken spectrum 

§r = A-^ forE.<Eb (4) 

, 1 , ^ ^ 
A— for < Eb 



dE^ El 



and a continuous spectrum = A^. We used the parametrization for A^^j^ given in 
Ref . which describes the response of the IceCube detector after all backgrounds have been 



rejected (this is achieved by quality cuts referred to as "level 2" cuts in Ref. [2J]) and we fix 

Irain = 300 m. 

The figure illustrates the challenge in experimentally determining the presence of an effect 
which depends on the original neutrino energy based on the observation of a muon track 
event. We can furthermore infer the even bigger challenge for the case of establishing (or 
constraining) the violation of Lorentz invariance from the observation of the arrival times of 
two muon events with different energy from a single GRB, a case we discuss next. 

IV. RESULTS AND CONCLUSIONS 

In principle one could test the hypothesis of LIV from the observed arrival times, ti and 
t2, of two neutrino events from a single GRB by confronting 

Aiobs = h — tl = '^tsource + ^t2,LIV ~ ^tl,LIV (S) 

with At given in Eq.([2]). Here we denote by Atgource the time difference of production of 
the two neutrinos at the source and not at arrival. It is not directly observed. Assuming we 
can infer Atgource, then for n = 1 

1 Fi(z)E,,2-i^.,i 

A J. A J. P 

-no ^Planck 



where we have defined 



F„(z) = / = . (7) 
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Log.o (EyCeV) Log,o(EyGeV) 



FIG. 1: Probability distribution for the parent neutrino energy of a muon track event which deposits 
an energy E^!:^ at the detector for different values of E'^'''^ = 10^ 10^ 10^, 10^, 10*^, and 10^ GeV 
as labeled in the figure. The left panel was obtained for a broken neutrino spectrum (Eq.([3])) with 
El, = 10^ GeV while the right panel corresponds to a continuous -E"^ spectrum. 

Numerically for Vt^ = 0.25 and = 0.75 0.28 < Fi{z) < 2.5 and 0.35 < F2{z) < 8.5 for 
0.5 < 2; < 5. 

However, as illustrated in FigiU inferring E^^2 — E^^i from the measured values of E'^f 
and of two muon tracks is highly non-trivial. For example, for n = 1 one can determine 
the relevant probability distribution of the parent neutrino energy difference AE^ as: 

f{/\E^Efl E^i) = I dE,, I dE,J{E,^E'^J)f{E,^2Ki)m,2 - E^,, - AE,) . (8) 

(Equivalently one can build the corresponding probability distribution for any combination 

We show in Figj2] the value of /(AE^|Ejf , for E'^'^f = E^^ = 10^ GeV, E^'^f = 

W GeV and E"^^"^^ = 10^ GeV and Ejf = 10^ GeV and eJI' = 10* GeV. From the 
figure we see that, because of the relatively broad spectrum of the parent neutrino energy 
distributions, one cannot rule out relatively large parent neutrino energy differences even if 
the two observed muon track events deposit the same energy at the detector. Conversely, 
only for very different deposited energies one can establish that the parent neutrino energies 
were different with high CL. 
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FIG. 2: Probability distribution for the parent neutrino energy difference of two muon track events 
which deposit an energy E'^^f and E'^^^ inside the detector for different values of E'^^f and E'^^^f as 
labeled in the figure. This distribution is for a E~'^ neutrino spectrum. In all cases f{A)dA = 1 
but we are showing /(A) as a function of log(A) hence the apparent difference in normalization of 
the curves. 

More quantitatively for a given pair of observed events with E'^^f and -E^^l" from a single 
GRB we can determine: 

f{A\E'^';f,E^J)dA, 

-oo 

POO 

Prob(E,,2 - ^.,1 > Amin) = / /(Al^Jf , dA . (9) 



From Eqs.(l6]) and ([9]) one can compute the probability for E,max > ^ > C,min which corresponds 

<:i E^l 



to the observation of a pair of events from a single GRB with E'^'^^ and -E'!^f and arrival 



times ti and t2- 

As illustration we show in Figsj3] and H] the ranges of ^ which can be inferred from the 



observation of two muon track events with E^T = 10^ GeV and ET^ = 10« GeV and with 
E^J = W GeV and = lO^ GeV respectively, ClS db function of the ratio of the arrival 
time difference and the time difference at production. 
From the figures we see that: 



pdep 



What one observes is in all cases compatible with no LIV if the model predicts 
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FIG. 3: Allowed ranges of the LIV scale factor ^ with probability bigger than 10~^ and 10^^ from 
the observation of two track events from a single GRB with E'^'^f = 10^ and -E^'f' ~ GeV as a 
function of the ratio of between the observed time difference and the time difference at production 
in the GRB. This is shown for an assumed continuous E^'^ neutrino spectrum. 

source = ^iobs, heiice the asymptotic peaks which reach ^ = oo and ^ = — oo in 
the upper and lower right panels at At source = Atobs- However, as can be seen in 
Fig. [3l if the deposit energy of the two events are not different enough, the allowed 
range of ^ include very low values. Thus no meaningful bound on the scale of LIV can 
be set. 

• When the deposited energies of the two events are very different and one assumes 
that Atobs = Atsource, ouc cau set a lower bound on the scale of LIV. For example for 
£!jP = 10^ GeV and = 10^ GeV we find at 99.9% probability that 

1^1 > 10«Fi(.)-^ (10) 

source 

• In the absence of an argument for Atsource = Atobs , the conclusion is that one has to 
be observing LIV. However, as seen in Fig. [3l if the two events deposit energies which 
are not different enough the observation is compatible with a very wide range of scales 
of LIV and with both signs of the effect. 
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FIG. 4: Same as FigElbut for E'^J = 10^ GeV and E'^^'^^ = 10^ GeV. 

When the deposited energies of the two events are very different the determination of 
^ becomes better as imphed by the narrow width of the allowed regions in Fig. HI Also 
one sign is favored: it is more likely that ^ has the opposite (same) sign than Atsource 

if the Atobs > Atsource {Atobs < Ats 



^source I 



Finally, the estimates assume negligible background from atmospheric neutrinos. One 
expects of the order of 200 muon track events a year with E^*^^ > 10^ GeV arriving from 
any direction with zenith angle cos^ < 0.2 [2,]. This means that, on average, the expected 
number of atmospheric background events in an angular bin of A^ioxio ^ 3 x 10~^ sr and 
in a time interval of Atobs is 

^atm,v^ r^dep y ^g4 Q^y^l ^ 3 X 10"^ (ll) 

^ min 

So as long as the observed time window is at most of the order of months the atmospheric 
background can be safely neglected even if one of the events has deposited energy as low as 
10^ GeV. 

In summary revealing Planck-scale physics by comparing the arrival times of muon track 
events from a single GRB at IceCube may be possible, but will be challenging. 

Observing the less frequent i^e and i^r events could present us with a solution to the 
energy measurement problem. Because of oscillations over cosmic distances the GRB flux 
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will be equally distributed over the three neutrino flavors. Neutrino telescopes detect the 
Cherenkov light radiated by secondary particle showers produced by neutrinos of all flavors. 
These include the electromagnetic and hadronic showers initiated by v^. and as well as 
by neutral current interactions of neutrinos of all flavors. Because the size of these showers, 
of order 10 m in ice, is small compared to the spacing of the PMTs, they represent, to a 
good approximation, a point source of Cherenkov photons radiated by the shower particles. 
These trigger the PMTs at the single photoelectron level over essentially spherical volume, 
slightly elongated in the direction of the initial neutrino, whose radius scales linearly with 
the shower energy [lo|. 

Whereas the smaller flrst-generation telescopes mostly exploit the large range of the 
muon to increase their effective area for i/^, kilometer-scale detectors can fully exploit the 
advantages associated with the detection of showers initiated by Vf, and Vt-: 

1. They are detected over both Northern and Southern hemispheres. 

2. The background of atmospheric neutrinos is signiflcantly reduced. At higher energies 
the muons from vr decay, the source of atmospheric i^e, no longer decay and relatively 
rare K-decays become the dominant source of background electron neutrinos. 

3. Vr are not absorbed, but degraded by energy in the earth. 

Critical in this context is that, depending on energy, electron and tau neutrinos deposit 
0.5-0.8% of their energy into an electromagnetic shower initiated by the leading flnal state 
lepton. The rest of the energy goes into the fragments of the target that produce a second 
subdominant shower. For ice, the Cherenkov light generated by shower particles spreads 
over a volume of radius 130 m at lOTeV and 460 m at lOEeV, i.e. the shower radius grows 
by just over 50 m per decade in energy. The measurement of the radius of the lightpool 
mapped by the lattice of PMTs determines the energy and turns neutrino telescopes into 
total absorption calorimeters. Note that even a contained "direct hit" by a 10 EeV neutrino 
will not saturate a km^ detector volume. So, even for EeV neutrinos, IceCube will not 
saturate and their energy spectrum can be measured. 

The drawback of the use of shower events are the lower statistics and the worse angular 
resolution. Although, in principle, GRBs represent a background-free signal in neutrino tele- 
scopes because of the clustering in time and in direction, the event rates are still uncertain. 
They will exceed 10 events per kilometer square year, in the case where the baryon loading 
in GRB flreballs accommodates the observed cosmic rays. With the low event rates antici- 
pated doing the Planck science with showers may not be a practical solution. Furthermore 
their worse pointing will also difficulty background rejection. 

In the end we will have to await the actual observations to draw any conclusions regarding 
the reality of observing Planck physics with neutrino telescopes. We have argued that the 
existing back-of-the-envelope estimates are unlikely to be realized and have described the 
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realistic background on which the possibihty of observing Planck physics will have to be 
pursued. The problem is undoubtedly important given the absence of clear alternatives. 
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